MATCHING CONTROL LAWS FOR A BALL AND BEAM SYSTEM 
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1. THE MATCHING CONDITION 
This note presents an application of the method 
developed by Auckly, et al. (2000), to stabiliza- 
tion of a ball and beam system. The results are 
fully described in (Andreev, et al., (2000), Auckly, 
Kapitanski (2000), and Auckly, et al. (2000)). An 
experimental comparison of a linear control law 
versus the nonlinear control laws described here 
will be given in the full paper, (Andreev, et al., 
(2000)). 

Let Q denote a configuration space. Let g £ 
T(T*Q®T*Q) be a metric. Let c, / : TQ -> TQ be 
fiber-preserving maps. We assume that c(— X) — 
—c(X). Let V : Q — ► R. The differential equation 
that we consider is 



V^7 + c(j) + gradyV = / (7). 



(1) 



Let P £ T(T*Q ® TQ) be a g-orthogonal projec- 
tion. We consider the situation where a constraint 
P(f) = is imposed. A system is called underac- 
tuated if P ^ 0. 

Several recent papers propose to find control in- 
puts so that the closed-loop system (Q) would 
have a natural candidate for a Lyapunov function 
(Bloch, et al. (1998), Hamberg (1999), and van 
der Shaft (1986)). Auckly, et al. (2000) intro- 
duced the following matching condition and char- 
acterization of matching in terms of linear partial 
differential equations. A control input, /, satisfies 
the matching condition if there are functions g, c, 
and V so that the closed loop equations take the 
form: 



V^7 + c( 7 ) + grad V = 0. 



(2) 



The motivation for this method is that H = 
59(7, 7) + ^(7) is a natural candidate for a Lya- 
punov function because dH/dt = —g(c(j),j). 
A straightforward computation shows that, the 
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matching condition is satisfied if and only if 

P{V X X - V X X) = 0, (3) 



(4) 



P{grad 1 V — grad^V) = 0, 
P(c(X)~c(X)) = 0. 

Equation (|3|) is a system of non-linear first order 
PDE's for g. It is perhaps surprising and pleas- 
ing that all of the solutions to (||), (||) may be ob- 
tained by first solving one first order linear system 
of PDE's and then solving a second set of linear 
PDE's. This is accomplished by introducing a new 
variable, A, by g(X, Y) = g(XX, Y). 

Theorem 1 The metric, g, satisfies (Q) if and 
only if A and g satisfy 
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(5) 



In the special case of a system with two de- 
grees of freedom, it is possible to write out the 
general solution to this set of differential equa- 
tions. Following Auckly, Kapitanski (2000), ex- 
press the underactuaded subspace as the span of 
a unit length vectorfield, FX. Choose coordinates 
x 1 , x 2 so that PX = and write XPX = 

cr^p- + Mgf^- F° r the A-equation, (|5|), to be con- 
sistent the following compatibility condition must 
hold: d([ll,2]/x)/&r 2 = d([U,2] rf/dx 1 . Start- 
ing with this equation and working backwards, all 
of the equations may be solved via the method of 
characteristics. 

2. THE BALL AND BEAM SYSTEM 




Fig.l. Nonlinear mechanical system. 

As an application of our method consider the sta- 
bilization problem for the ball and beam system 



described schematically in figure 1. One can ex- 
press a as an explicit function of 9. After rescal- 
ing, the kinetic energy of the system is given by: 

T = \s 2 + a'sO + \{a i+ L + (af) 2 

and V = a5sin(0) + (s + a6)sin(a), where the 
afe are dimentionless parameters. The projection, 
P = (ds + a'dO) ® d/ds, so the control input u is 
related to / in (1) by / = (ud9)K The resulting 
equations of motion are 

s + a'6 + (a" - |sc/ 2 )<9 2 + sin(a) = 

a's + [a 4 + (a 3 + §s 2 )a' 2 ]<9 + ha' 2 ss9 
+ (a 3 + \s 2 )a'a"9 2 + a 5 cos 9 

+(ae + s) cos(a) a' + aj9 = u , 

where corresponds to inherent dissipation. 

The general solution to the matching equations is 

g ll {s,6)=^{a) (%M)) + 10 f %- ) 

Jo ViivW 'VPS 

9i2 = -(.9n - vgu), 922 = -{gu - 0-512), 
V(s,9) = w(y) + 5 (y + s ) J" d(p 

where y = ip(a)s — sq + J Q Q ip( T ) dr, ip(a) = 

/11(a) — and /ii, h, and w are arbitrary 

functions. Also, c 1 = — a'c 2 , where c 2 (s, 6, s, 9) is 
an arbitrary function which is odd in s and 9. The 
final nonlinear control law is u = u g + uy + u c , 
where u g = g{V^ ~ V^7,^), u v = % - 
g(grad 7 V, §g), and u c = a 7 8 - g(c(j), ■§$). Using 
H as a Lyapunov function, we obtain the following 
conditions that guarantee local asymptotic stabil- 
ity of the equilibrium: det(<?(0)) > 0, tr(<?(0)) > 0, 
det(5c(0)) > 0, tr(gc(0)) > 0, det(D 2 l>(0)) > 0, 
and tr(L> 2 \>(0)) > 0, 

Another way to check local asymptotic stability 
is to find the poles of the linearized closed-loop 
system. It is a theorem (Andreev, et at (2000), 
Auckly, Kapitanski (2000)) that any linear full 
state feedback control law can be obtained as a 
linearization of some control law in our family. 

A good stabilizing control law will produce a large 
basin of attraction, send solutions to the equilib- 
rium in a short period of time, and will require 
little control effort. It is, unfortunately, not clear 
how to quantify these goals. 

We have done some numerical simulation of vari- 
ous control laws in our family. We always pick the 



arbitrary functions in our nonlinear control law 
in such a way that the linearization at the desired 
equilibrium, u Mn = a s +K bp (s-s )+K ap 9+K bd -s+ 
K a d9, is exactly the linear state feedback control 
law provided by the manufacturer of a commer- 
cially available system (Apkarian, (1994)). The 
numerical and experimental response of the sys- 
tem to various initial conditions will be recorded 
in the full version of the paper. 

3. CONCLUSION 
We believe that nonlinear control laws have the 
potential to achieve better performance than lin- 
ear control laws. There are, however, several sub- 
tle questions which must be resolved before non- 
linear control laws may be fully exploited in prac- 
tice. The first question is how to quantify per- 
formance. The second question is how to pick a 
control law which will come close to optimizing 
performance. One interesting idea is to restrict 
attention to a class of control laws which generate 
a closed loop system of a special form. The hope 
is then that it will be easier to quantify the per- 
formance of such systems. We have shown that, 
in many situations it is possible to find all control 
laws which will result in a closed loop system of 
the form (^). 
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